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Abstract 

The pion re-scattering operator for pion production, derived recently in time-ordered perturba- 
tion theory, is compared with the one following from the simple S-matrix construction. We show 
that this construction is equivalent to the on-shell approximation introduced in previous papers. 
For a realistic A^A^ interaction, the S-matrix approach, and its simplified fixed threshold-kinematics 
version, work well near threshold. 
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I. INTRODUCTION 



The detailed analysis of the irreducible pion re-scattering operator was recently per- 
formed for the reaction pp ppn^. The pion re-scattering is certainly part of the pion 
production mechanism, but its importance relative to other contributions varies consider- 
ably dependent on the approximations made in evaluation of the effective operators (see 
Ref. P] and references therein). The nature and extent of this uncertainty are re-examined 
in this paper. 

To this end we deal with retardation effects in the exchanged pion propagator, i.e., its 
energy dependence, as well as with the energy dependence of the vrN scattering amplitude in 
the vertex, from which the produced pion is emitted. While the approximations employed in 
a previous paper llj yield rather different results, we show here that the deviations between 
them are significantly reduced when the approximations are applied consistently in the whole 
effective operator. 

We also show that the S-matrix approach, which has been successfully used below pion 
production threshold, yields also above threshold results rather close to those obtained with 
the energy-dependent operator following from time-ordered perturbation theory. 

The paper is organized as follows: after this Introduction, section II describes the S- 
matrix technique for deriving effective nuclear quantum-mechanical operators, section III 
describes the results and section IV presents a summary and conclusions. 



II. PION RE-SCATTERING OPERATOR 



To derive the effective pion production operators, and other effective nuclear operators in 
general, one starts from the relativistic (effective) Lagrangian written in terms of hadronic 
fields. The interactions mediated by meson exchanges before and after the production 
reaction takes place are included in the effective NN (and nucleon-meson) interaction, while 
from the irreducible parts connected to the reaction mechanism (e.g., pion production) one 
obtains effective operators, whose expectation values are to be evaluated between the initial 
and final nucleonic wave functions. One aims is to get such effective operators consistent 
with the realistic description of the A^A^ interaction, which can then be used in studies 
of the corresponding reactions not only on the simplest (one or two-nucleon) systems, but 
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preferably also on heavier nuclei. 

The covariant techniques based on the Bethe-Salpeter equation or its quasipotential re- 
arrangements are these days practically manageable only below meson production threshold. 
However, above the threshold the dressings of the single hadron propagators and interaction 
vertices via the meson loops have to be included explicitly. For this reason the construction 
of the production operator is so far realized mostly in the Hamiltonian quantum-mechanical 
framework (usually non-relativistic, or with leading relativistic effects included perturba- 
tively within the decomposition in powers of p/m, where p is typical hadronic momentum 
and m is the nucleon mass). 

The derivation of the nuclear effective operators below the meson production threshold 
within the Hamiltonian framework - leading to hermitian and energy independent A^A^ and 
3A^ potentials and conserved e.m. and partially conserved weak current operators - can be 
done in many different ways (see discussion in Ref. P] and references therein). At the non- 
relativistic order the results are determined uniquely. As for the leading order relativistic 
contributions, they were shown to be unitarily equivalent. The unitary freedom allows to 
choose the A^A^ potentials to be static (in the cm. frame of two-nucleon system) and identify 
them with the successful static semi-phenomenological potentials. 

Also above the threshold the static limit is commonly employed, since more elaborate 
descriptions which include the mesonic retardation and loop effects are technically consider- 
ably more complex 0, |^ , especially for systems of more than two nucleons. Both the static 
approaches and the ones including "retardation" typically consider contributions of several 
one-meson exchanges and the potentials are fitted to describe the data. It is therefore diffi- 
cult to assess how well do they approximate the covariant amplitudes (corresponding to the 
same values of physical masses and coupling constants) which are so far outside the scope 
of existing calculations schemes, but which we believe do provide in principle the consistent 
description of the considered reactions. 

Thus, the ultimately exact approach to the description of the pion production (and 
in particular of the pion re-scattering contribution) would be either the covariant Bethe- 
Salpeter or quasipotential frameworks (extended above the pion threshold) or the quantum 
mechanical coupled-channel technique including retardation. In these approaches one has 
to treat the non-hermitian energy-dependent A^A^ interaction (fitted to the data also above 
pion production threshold) and consider the effects of renormalization of vertices, masses 
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FIG. 1: Feynman diagrams for pion re-scattering. The pion field is represented by a dashed hne, 
the NN interaction by sohd double line and the nucleons by solid lines. 

and wave functions via meson loops. 

In this paper we rather (following Refs. numerically estimate the range of the pre- 
dictions from several commonly used simplifying approximations, and compare them to the 
result obtained from the reduction of the corresponding covariant Feynman diagrams for 
the pion re-scattering operator. This reduction coincides with the time-ordered perturba- 
tion theory 



A. Factorization of the effective re-scattering operator 

In a previous paper [1] we made the connection to the usual representation of the pion 
re-scattering operator for non-relativistic calculations by following the approach of Refs. P| . 
We started from the covariant two-nucleon Feynman amplitudes including final and initial 
state interaction (FSI and ISI, respectively), shown in Figs. ^ and^. 

To obtain the effective re-scattering operator the negative energy contributions in the 
nucleon propagators (to be included in the complete calculations) were neglected. By inte- 
grating subsequently over the energy of the exchanged pion the resulting Feynman ampli- 
tudes were transformed into those following from the time-ordered perturbation theory. We 
have shown in Ref. flj that the irreducible "stretched box diagrams" (i.e., those with more 
than one meson in flight in the intermediate states) give a very small contribution and can 
be therefore also neglected. Thus, the full covariant amplitude is in the lowest order Born 
approximation well approximated by the product of the NN potential and the effective pion 



re-scattering operator, which can be extracted from these time-ordered diagrams (Fig. |2I). 

The effective pion re-scattering operator was in Ref. |l| factorized into an effective pion re- 
scattering vertex / and an effective pion propagator G^. For the diagram with FSI (Fig. |21) 
this factorization reads 

MfsI = / 7^ -ET—B Ors , (1) 

J {2Tiy Fi + F2 - coi - UJ2 + le 



E2-UJ2- t^TT El - uji- - 10^ 



\~fG^. (2) 



/ = — [(^1 -uJi-E^- uj^)fM + {E2 -UJ2- uj^)f{-uj^)] , (3) 

(El -Ui- E^- U^) {E2 -UJ2- LO^) ' 



1 1 

+ 



UO2 —ijJa Fi —UJl — UJo 



(5) 



where, adopting the notation of Ref. g' is the momentum of the exchanged pion, ul = 
ml + q''^ is its on-mass-shell energy, /(cUtt) is the product of the vrA^ amplitude with the nNN 
vertex (as in Ref. jl| the standard xPT re-scattering vertex is employed here), Ei,uJi,Fi are 
the on-mass-shell energies of the i-th nucleon in the initial, intermediate and final state, 
respectively, E.„ is the energy of the produced pion, Ei + E2 = Fi + F2 + E.,,. 

The inclusion of some pieces of the integrand of Eq. ((U) into the propagator and 
of others in the modified vertex / is somewhat arbitrary. The appearance of the unusual 
effective propagator G^^ and the effective vertex / is the result of combining two time-ordered 
diagrams with different energy dependence into a single effective operator. 

The A^A^ interaction is in Fig.|2land in Eq. (0) simulated by a simple a-exchange potential. 
Though not realistic, this interaction suffices for model studies of approximations employed 
in derivations of the effective pion re-scattering operator, as done in references Q]. Since 
some results do depend on the behavior of the A^A^ scattering wave function, in particular 
in the region of higher relative momenta, we perform our calculations (as in Ref. ^|) also 
with V(j replaced by a full A^A^ T-matrix, generated from realistic Bonn B potential. 

We note that the meson poles are not neglected in the integration over the energy Qq 
of the exchanged pion, which generates Eq. (^. A result similar to Eq. (0) can also be 
obtained for the amplitude with the initial state interaction (ISI). The two amplitudes differ 
however in the contribution from the pion poles to the remaining integration over the three- 
momentum. For the amplitude with FSI there are no such poles. However, for the ISI 
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FIG. 2: The two time-ordered diagrams for FSI considered here, additional stretched box diagrams 
are neglected. 

case there are values of the exchanged pion three-momentum for which the propagator 
has poles. These poles have been considered in all our numerical calculations for the cross 
section. As we will see, they are one of the main reasons for deviations between several 
approximations and the reference results calculated from Eqs. (fTJEl). 

It is worth mentioning that although the FSI and ISI diagrams graphically separate the 
A^A^ interaction and the pion re-scattering part (when the stretched boxes are neglected), 
they do not define a single effective operator (as a function of nucleon three-momenta and 
the energy of emitted pion). Since in these time-ordered diagrams energy is not conserved 
at individual vertices, each of these diagrams defines a different off-energy shell extension 
of the pion re-scattering amplitude. This is an unpleasant feature, since one would have to 
make an analogous construction for diagrams with both FSI and ISI. Moreover, one would 
have to repeat the whole analysis for systems of more than two nucleons. Only after the 
on-shell approximation is made (in the next subsection), the pion re-scattering parts of FSI 
and ISI diagrams coincide and one can identify them with a single effective re-scattering 
operator. 



B. Re-scattering operator in the S-matrix technique 

The S-matrix technique is a simple prescription to derive the effective nuclear operators 
from the corresponding covariant Feynman diagrams j^. For electromagnetic operators 
and also for A^A^ and 3A^ potentials the S-matrix approach reproduces the results of more 
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laborious constructions, based on time-ordered or non-relativistic diagram techniques. 

The two-nucleon effective operators are by definition identified with the diagrams de- 
scribing the irreducible mechanism of the corresponding reaction. The only exception are 
the nucleon Born diagrams from which the iteration of the one-nucleon operator has to be 
subtracted. The operators of the nuclear e.m. and weak currents, as well as the pion ab- 
sorption operators and nuclear potentials, are obtained by a straightforward non-relativistic 
reduction of the corresponding Feynman diagrams, in which the intermediate particles are 
off-mass-shell and energy is conserved at each vertex: therefore the derived effective oper- 
ators are also defined on-energy-shell. The nuclear currents and other transition operators 
are defined to be consistent with a hermitian energy independent A^A^ potential, which has 
the usual one boson exchange form employed in realistic models of A^A^ interaction, and 
can be used also in systems of more than two nucleons. This approach is well defined and 
understood below the meson production threshold, but as a simple tool it is employed also 
above the threshold, for instance in Refs. 3, namely to derive the Z-diagram operators. 

For the a exchange potential the S-matrix technique in the lowest order of non-relativistic 
reduction yields 

y yon ^ \ /g\ 

with A = Ai = — A2 and Aj = e'- — e^; and being the on-shell energies of the i-th nucleon 
after and before the meson exchange, respectively. As pointed out above, this defines the 
potential only on-energy-shell. However, the Lippmann-Schwinger equation and even the 
first order Born approximation to the wave function require the potential off-energy-shell. 
The extended S-matrix approach defines the most general off-energy-shell continuation of 
Va- as a class of unitarily equivalent potentials parameterized by the "retardation parameter" 
u. The particular choice u = 1/2 leads to the static potential in the A^A^ cm. frame. This 
choice corresponds to the substitution A^ = (Ai — A2)^/4. Most realistic A^A^ potentials, 
namely those fitted to the data below pion threshold, in particular the Bonn B potential 
used in this paper, are energy-independent and static in the nucleon cm. frame and can be 
therefore considered to be consistent with this construction. 

For the pion re-scattering diagram S-matrix prescription leads to a single effective oper- 
ator (both for FSI and ISI diagrams) of the form: 



where = e'2 — e2 = ei — e'l- 

Let us finally introduce approximations to the time-ordered PT result given by Eqs. 
and explain how the S-matrix technique fits in. The first of them is the on-shell approxi- 
mation, in which the two nucleons in the intermediate state are put on-energy shell. That 
is, we put Ui + UJ2 = Fi + F2, which also implies Ei + E2 = ui + UJ2 + E^,. For the scalar 
potential in the FSI diagram this leads to ©, where now A = Fi — ui = uj2 — F2 is the 
energy transfer in the corresponding vertices. 

For the re-scattering operator (0) the on-shell replacement implies 



rs 



(8) 



E2 — LO2 — lOt^ UJ2 — E2 — o;. 

Clearly, the S-matrix pion re-scattering operator O^^ follows from 0°.l if one assumes the 
energy conservation at each vertex. The exchanged pion is then no longer on-mass-shell and 
we have to replace /(tUvr) — > /(-E'2 —^2) and /(— cj^) /(-E'2 — ^2)'- in the first time-ordered 
diagram the virtual pion is entering the re-scattering vertex and in the second one it is 
emitted from this vertex (as defined on Fig. El). 

The on-shell approximation as introduced in Refs. 0, 0] actually coincides with the S- 
matrix approximation defined above. The re-scattering operator ((7j) can be obtained directly 
from and (ji} by the substitutions following from the on-energy-shell prescription and 
the energy conservation in individual vertices uOt^ = E2 — 002 = — {Ei — cui — E.,^) as follows: 
A 1 -2 {E2 - 0^2) X / {E2 - ^2) + X / (^2 - E2) ^ 

2(^2-^2) {E2 - 002-00^) {002 -E2-UJ^) 

In Ref. the extra kinematical factors in (jS)) multiplying the function / (tu^) were inter- 
preted as form factors and kept unaltered, i.e., the substitution above was made only in 
and f{uJn) of Q, not in the kinematical factors included in the function /. 

In equation the effective pion propagator Gt^ is seen to take its Klein-Gordon form: 

G.^GT = iniE2-co2r-col]. (10) 

The replacement (jiup does not significantly alter the results, as we will show below, but the 
corresponding substitution alone in the pion re-scattering vertex /(icj^) f{±{E2 — coioV) 
present in Eq. (jH} leads to a large enhancement of the amplitude (P). According to Ref. p] 
it increases the cross section by almost factor of 3. Since the splitting of Ors is not defined 
unambiguously, in the work reported here the on-shell replacement is made in the whole 
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re-scattering operator. We demonstrate in this paper that making the on-shell replacement 
in the whole operator leads to a significantly smaller deviation from the reference result, 
compared to the replacement /(ic^Tr) f{±{E2 — 002)) involving / only . 

In our previous paper we considered also other approximations (besides the on-shell 
one): the so-called static and fixed threshold-kinematics approximations, defined by the 
replacement of the energy of the exchanged pion E2 — 002 by zero and 171^^/2 (its threshold 
value), respectively. The static approximation was in Ref. |l| considered only for Gt^, the 
fixed threshold-kinematics one for and also for f{±ujjj) (the additional kinematical factors 
in / were again kept unchanged). In this paper we again make these approximations in the 
full re-scattering operator ((Tj), replacing f2 ^ and fl m.„/2, respectively. For Vo- the 
static approximation is defined by A ^ 0. 



III. RESULTS 



For numerical calculations we consider the A^A^ —>■ {NN)7t transition in partial waves 
^Pq (^So)so. Amplitudes and cross sections are evaluated both with the simple interaction 
Va and with the Bonn B potential. We test the S-matrix prescription for the re-scattering 
operator (|7j) and also the fixed threshold-kinematics and the static approximations discussed 
in the previous section. Besides, in order to compare to the previous papers we include also 
the results for the on-shell fixed threshold-kinematics and the static approximations 
for the effective pion propagator G^. 

In Fig. El we show that the amplitudes with the S-matrix operator (dotted line with 
crosses on the upper panels) are the closest to the reference result (solid line). Using the 
same approach both for the operator and for V„ increases slightly the gap from the reference 
result (dashed-dotted line versus solid line on the upper left panel). The fixed threshold- 
kinematics version of Ors, denoted as O^'', works well for small values of the excess energy 
Q = 2E — 2M — Et^, but starts to deviate rapidly with increasing Q (dotted line in the upper 
panels of Fig.E}. The static approximation for the re-scattering operator (O**) overestimates 
significantly the amplitude ((H) (dashed versus solid lines in the upper panels). 

From the lower panels on Fig. |H1 one sees that all considered approximations taken only 
for the effectivepion propagator do not differ much from each other. This had already been 
found on Ref. It means that the choices for the energy of the exchanged pion in the 
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effective propagator alone are not very decisive (solid line versus dotted, dashed and cross- 
dotted lines). We notice however that there is a considerable deviation (dependent on the 
A^A^ interaction employed) of all these approximations from the reference result. 

In order to understand this we considered the expansion of the effective pion propagator 
in Eq. (@)) in terms of an "off-mass-shell" dimensionless parameter y: 

y = — ) (11) 

UJl- UJ2 + -C/vr 

which measures the deviation of the total energy from the energy of the intermediate state 
with all three particles on-mass-shell. This Taylor series expansion gives insight on the small 
effect of retardation effects in the propagator, and it reads 



, -2E + + UJl + U02) 

1 + — 72 + 



ETr+U)l-U)2 \ _ , ,2 
2 ) ^-n 



(12) 



^Tay 



where G^ly has the form of the usual Klein-Gordon propagator. 

We notice here that in the case of the ISI amplitude, the representation of the pion 
propagator by its Taylor series, the first term of which is G^ly, fails due to the presence 
of a pole in the propagator. 

Figure m compares the first four terms Gjl^y{i = 1, ...,4) of this expansion with the full 
effective propagator in Eq. as a function of the two nucleon relative momentum q^, for 
two different values of the excess energy Q = 2E — 2M — m^. The convergence of the series 
demands at least 4 terms. Besides, as expected, this convergence is momentum-dependent. 
We have also compared the first term of this expansion with the already considered on-shell, 
fixed threshold-kinematics and static approximations for the pion propagator. These results 
are shown on Fig. El We realize that all these approximations are very near to the 1st order 
term of the Taylor series. The corrections arising from higher order terms in the expansion 
are negligible only for low momentum transfer, more precisely in the range < lOOMeV. 

The deviations of G**, G^'' and from the effective propagator G^r given by Eq. (j3} 
cannot explain the relatively large deviations obtained on the bottom-left panel of Fig. |21 
between considered approximations and the reference result. These deviations follow from 
the ISI contribution. The weight of the ISI term depends on the NN interaction employed. 
It is comparable to the FSI term for (for which the deviations are large, as seen on the 
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bottom-left panel of Fig. EI), but it is much less important for the full Bonn B potential (and 
therefore the corresponding deviations on the bottom-right panel of Fig. El are indeed much 
smaller) . 

All the findings for the amplitudes manifest themselves also in the results for the cross 
section. We show in Fig. IHlthe effects of the considered approximations on the cross section, 
first taking only the FSI contribution. On the left panel the amplitude includes for the 
A^A^ interaction, on the right panel the Bonn B potential is used. The curves compare the 
reference result (solid line in all panels) with the S-matrix results (upper panels) and their 
fixed threshold-kinematics version (lower panels). The S-matrix approach (dashed line) is 
the closest to the reference result (upper panels of Fig. 

For the case of the A^A^ interaction described by V„ we also show the result following 
from the S-matrix prescription applied to the A^A^ interaction (dotted line on left panels in 
Fig. in)). For the fixed threshold-kinematics versions (bottom- left panel) the deviations from 
the reference result increase more pronouncedly with the excess energy Q, as expected. The 
approximations for the energy of the exchanged pion taken in the pion propagator and 
in the re-scattering vertex f{uJn), but not in the kinematical factors of (jH)), overestimate the 
cross section by a factor of 5 (solid line with bullets). 

Finally, we present on Fig. [7| the comparison between the approximated total cross sec- 
tions with both FSI and ISI included. The approximation dictated by the S-matrix approach 
(dashed and dotted lines on the upper panels of Fig. Ej) is clearly seen as the best one. For 
the Bonn potential calculation, it practically coincides with the reference result. As shown 
in the previous section, this procedure amounts to extend the on-shell approximation, used 
in Ref. jl| for and /(t^vr) alone, also to the multiplicative kinematical factors showing up 
in the operator / (see Eqs. ((21111) )• 

To conclude we notice, moreover, that for the realistic A^A^ interaction, the difference be- 
tween the S-matrix approach (upper right panel on Fig. |2|) and its fixed threshold- kinematics 
version (lower right panel of the same figure) is not very important near threshold, provided 
the excess energy does not exceed ~ 30 MeV (Q/m^ ~ 0.2). 
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IV. CONCLUSIONS 

1) The usual approximations to the effective pion propagator [Jj, |6|, |9|, obtained from a 
quantum-mechanical reduction of the Feynman diagram describing the re-scattering, which 
are rather close to the first order term of a Taylor series in a parameter measuring off-mass- 
shell effects in the intermediate states. The series converges rapidly for the FSI amplitude 
near threshold. As a consequence, retardation effects are not decisive in the pion re-scattering 
mechanism near the threshold energy for pion production. 

2) As for the pion energy in the vriV re-scattering amplitude, the on-shell approach when 
used only in /(cUtt) overestimates significantly the reference result. Nevertheless, and this 
is the key point of this paper, this deviation is dramatically reduced if the approximation 
coming from the S-matrix approach is used consistently in the whole effective operator. 
This procedure amounts to extend the on-shell approximation used in Ref. for and 
/(cij^) to the full operator /, including kinematical factors which differently weight the two 
dominant time-ordered diagrams. The amplitudes and cross sections obtained with the S- 
matrix effective operator are very close to those obtained with the time-ordered one in the 
considered kinematical region. 

The re-scattering operator of this paper for the neutral pion production in the isoscalar 
ttN channel thus indeed seems to be relatively unimportant: its enhancement reported in 
previous papers followed from inconsistent or too crude (static or fixed threshold-kinematics) 
treatment of the energy dependence of the effective operator. Our findings explain why the 
calculation of Ref. 6|, where the on-shell approximation is used, artificially enhances the 
contribution of the isoscalar re-scattering term. On the _otlier hand, importantly and in 

for the different production 



retrospect, our results support the choice done in Refs. [3,0, ^ 
operators considered. 

The re-scattering mechanism is filtered differently by other spin/isospin channels in pion 
production reactions. For charged pion production reactions the general irreducible re- 
scattering operator comprises also the dominant isovector Weinberg- Tomozawa term of the 
ttN amplitude, and its importance is therefore enhanced. Investigation of these channels 
within the approach of this paper is in progress. 
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FIG. 3: Absolute values of the FSI + ISI amplitude as a function of the excess energy Q = 
2E — 2M — Ej^ (in units of 171.,^). The right(left) panels correspond to the amplitudes with a 
exchange (Bonn B) for the A'^A^ interaction. The amplitudes are taken at the maximum pion 
momentum g^"^, determined by Q. The upper panels correspond to approximations for the whole 
operator O^si the lower panels to approximations for the pion propagator G,r only. The solid line 
(O^^) is the reference calculation. The dashed, dotted and cross-dotted lines correspond to the 
static, fixed threshold-kinematics and on-shell approximations, respectively. The corresponding 
operators are O^^, and G^^, G°^. All amplitudes were normalized by a factor defined 
by the maximum value of the reference result. 
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FIG. 4: Convergence of the Taylor expansion for the pion propagator G-,^ in the FSI amphtude, as 
a function of the two nucleon relative momentum (Eq. (|12() ). Left panel: at threshold; right panel: 
above threshold at maximum pion momentum for an excess energy Q of 0.5m,r- Bottom panels 
zoom into the region of low relative momentum. 
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FIG. 5: Approximations for G-^ in the FSI diagram and the first term of the Taylor series. Left 

and right panels with the same meaning as on Fig. ^ 

15 




FIG. 6: Effects of the approximations for the re-scattering operator O^s and for the effective pion 
propagator G-,^ as a function of the excess energy Q. The cross section curves shown correspond 
to the FSI amplitude alone. The upper panels correspond to the re-scattering operator O given 
by Eq. Q and the lower panels to fixed threshold-kinematics approximation. The solid line is 
the reference calculation The dashed line is the S-matrix calculation for the re-scattering 
operator O given by Eq. (jT)) (upper panels) and the fixed threshold-kinematics approximation for 
(O (lower panels). The dotted line corresponds to take the S-matrix approximation not only for 
O, but also for the cj-exchange interaction The dashed-doted line corresponds to the on-shell 
(upper panels) and fixed threshold-kinematics (lower panels) prescriptions only for Gt^. The solid 
lines with bullets refer to the energy prescriptions taken for the propagator G,r and for /(wtt) in 
Eq. @, as in 1], but not for extra kinematical factors in /. All the cross sections were normalized 
with a factor defined by the maximum value of the reference result. 
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FIG. 7: The same of Fig. IHl but for the total (FSI+ISI) cross-section and considering only the 
approximations for O and G-,^. 
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